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We study a holographic model of a relativistic quantum system with a global U(l) sym- 
metry, at non-zero temperature and density. When the temperature falls below a critical 
value, we find a second-order superfluid phase transition with mean-field critical exponents. 
In the symmetry-broken phase, we determine the speed of second sound as a function of 
temperature. As the velocity of the superfluid component relative to the normal component 
increases, the superfluid transition goes through a tricritical point and becomes first-order. 

Qh! I. INTRODUCTION 

Recently, gauge/gravity duality has been used to model strongly-interacting systems in 

terms of a gravity dual. The most important case is the strongly-interacting quark gluon plasma 
created at RHIC. While a systematic treatment is still lacking, the Jv = 4 super- Yang- Mills plasma 
0^ | has been used instead, with some success. Most notably, the viscosity/entropy density ratio, 

universal among all theories with a gravity dual UL la, m , seems to be not too far away from the 
g; value extracted from experimental data flfl. 

The long-distance behavior of a black hole horizon is captured by the same hydrodynamic 
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equation governing the evolution of the plasma. As a result, from the gravity equations one can 
reconstruct the hydrodynamic equations, including transport coefficients. One can go beyond lead- 
ing order in the gravity equations and capture second-order corrections, in the process discoverin 



extra terms typically ignored in almost all implementations of the Israel-Stewart formalism [9. 10]. 
The approach has been extended to R-charged black holes, dual to fluids with chemical poten- 
tials QQ. 

In this work, we use black holes to study relativistic superfluids. Our work is a continuation 
of ref. [3], where a holographic model of a superfluid was constructed. Here, we investigate 
the behavior of the superfluid when a superfluid current flows through the system. It is known 
that in nonrelativistic superfluids there is a critical superfluid velocity, above which the superfluid 
phase does not exist. We discover exactly the same phenomenon in the relativistic superfluids at 
low temperatures: there is a first-order phase transition between the superfluid and the normal 
phase as one changes the superfluid current. The phase transition becomes second order at higher 
temperatures. 

We will also be discussing the superfluid hydrodynamics of such systems, in particular, the speed 
of second sound. A short discussion of hydrodynamics can be found in [141 ] . Hydrodynamics of a 
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system with spontaneously broken symmetries is different from hydrodynamics of normal liquids 
because the system contains long-range modes (Goldstone bosons) which must be included in the 
hydrodynamic equations. 



II. THERMODYNAMICS 



Let us first review the thermodynamics of quantum systems with a spontaneous U{1) symmetry 
breaking. 

For simplicity consider a system with one global U(l) symmetry, such as complex (ft 4 theory. 
We define the free energy F = — T In Z where 

Z = tTe~ pH , (1) 

and H is the Hamiltonian. While we eventually take the thermodynamic limit V — > oo, for now it 
is more convenient to think of the system as living inside a d — 1 dimensional torus T d_1 where 
the perimeters of the d—1 circles are Lj. 

When the U(l) symmetry is not spontaneously broken, the free energy F depends only on the 
temperature, and, if we do not take the thermodynamic limit, the Lf. F = F(T,Li). When the 
symmetry is spontaneously broken, the system contains a massless (pseudo) scalar field (p, which 
is compact: <p ~ (p + 2ir. In the example of the complex (ft 4 theory, ip is just the phase of (ft. If one 
thinks about <p as the elementary field, periodicity in ip space means that we are free to impose 
boundary conditions such that tp\ Xi =Li = <p\xi=o + 27m j. Thus there is a set of partition functions 
Z[n], where we integrate over all fields with boundary conditions specified by n. The usual ground 
state corresponds to rij = 0. 

As we shall see later, Z[n] does not exist in the strict thermodynamic sense. Physically, if one 
sets n-i / 0, the system will relax over time to a state with no winding. However, the time scale 
for this process may be long, and the winding state is practically stable. We assume stability. 

Define qi = di<p, an analogue of the topological charge density whose integral can distinguish 
various boundary conditions: 

/•Li 

q • dx = / qidx 1 = <p\ x% =u ~ ¥?k=0 = , (2) 
d Jo 

where Cj is the i'th spatial circle of the T d_1 . We can therefore impose the boundary conditions 
on (p by introducing Lagrange multipliers Aj which pick out the required values of nf. 

Z[n]=J D(fte' SE ^ = jDcftd\e- SEm exp^^A^^ q ■ dx - 2im^j . (3) 

We are taking advantage of the fact that the Euclidean path integral on a compact time circle of 
perimeter [3 = 1/T yields the thermal partition function at a temperature T. Having compactified 
the time direction, it is natural to ask what boundary conditions should be applied for tp in the 
time like direction. The answer to this question involves introducing a chemical potential, and we 
will return to it shortly. 
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We can write the field tp as tp{x) = onx % + 8tp{x) where on is constant, and 5<p(x) is periodic. The 
delta function in the path integral for Z[n) will pick out on = 2imi/Li, so knowing an is equivalent 
to knowing the boundary conditions for tp. In the thermodynamic limit, F depends on the Lj only 
through an overall multiplicative volume factor V and a» becomes a continuous variable which 
transforms as a vector under rotations. At zero temperature, the energy (for given n) is always 
minimized for tp = anx 1 because that's when the (Vtp) 2 term is smallest. So at zero temperature, 
the ground state can be uniquely characterized by the on- We assume now that the thermal state 
can be uniquely characterized by Oj and T, and define the partition function as 



Note that a is the equilbrium value of (V</?). In nonrelativistic superfluids 'Vtp is proportional to 
the superfluid velocity. We shall use the same terminology and call V</2 the superfluid velocity. Z[a] 
describes a stationary state with constant, non-zero superfluid velocity. The partition function Z 
is a scalar, and therefore can only depend on a 2 in the thermodynamic limit. 

Note that the state described by Z[a] is in fact metastable, rather than stable. This is because 
tp, as the phase of a complex field cp, is ill-defined when = 0; hence the winding numbers rij 
are not topological and can change with time. The process of "unwinding" can be visualized as 
follows. Imagine a state with only one of the rii equal to 1, and other winding numbers equal 
0. The state can be thought of as containing a (d—2) dimensional domain wall, across which tp 
changes by 2ir. However, in d > 2 dimensions, there is no topological conservation law which would 
ensure the stability of the domain wall, and the domain wall will decay through hole nucleation. 
In terms of superfluid hydrodynamics in 3+1 dimensions, nucleating a hole in the two-dimensional 
domain wall corresponds to producing a vortex loop. Production of vortex loops violates Landau's 
criterion for superfluidity [15(, and renders the state with non-zero a. metastable. 

In a theory with large N, like the one we will be considering, the rate of vortex loop production is 
exponentially suppressed, as both the energy of a vortex loop, and the action of the configuration 
describing vortex nucleations, are proportional to N. Therefore, provided that the volume and 
time scales are not too large, one can treat states with a nonzero superfluid velocity as thermal 
equilibrium states. 

We now would like to introduce a chemical potential fi and work in the grand canonical ensemble. 
We define the potential function f2 = F — fiQ = —PV = —TlnZ where 



and Q is the conserved charge corresponding to the U(l) symmetry. To introduce fi at the level 
of the path integral, it is useful to gauge the U(l) symmetry by coupling the system to a non- 
dynamical gauge field A^: 



We are temporarily ignoring the dependence of Z on the a. The one point function of the U(l) 





(5) 




(6) 
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current is then generated by W[A|: 

8W[A] 



<■/"(*)) . (7) 



In a system where the external field strength is zero F^ u = 0, we may pick a gauge in which is 
constant and it makes sense, given (|7|), to interpret Aq = [i as the chemical potential. 

Gauge transformations which send An —* A^ + d^X have a nontrivial effect on the compact 
scalar: if — > 95 + A. In the presence of the background gauge field, the pressure P can depend only 
on gauge invariant quantities such as D^ip = d^ip — A^. The grand canonical partition function 
thus becomes 

Z[A,a}= jDc/)dXe- SE ^ exp^^A^ (V<p - A- a) • dxj . (8) 



Like Z[a], by rotational symmetry 2 [a] can only depend on at in the thermodynamic limit. We are 



assuming that the pressure P = TV 1 ln2 depends on three parameters, P = P(T, fi, i(Z)j<£>) 2 ). 



We can now choose a gauge in which ip=0, so that {Dup) 2 = A 2 . The chemical potential n can 
be interpreted as the zero component of the gauge field, and we have P = P(T, Aq, Af). At zero 
temperature, the system is Lorentz invariant, provided one treats Aq and A^ as spurion fields that 
transform as a four-vector under the Lorentz group. As a result, the pressure can only depend 
on A 2 = A 2 —A\. At finite temperature, the pressure depends on Aq and A 2 separately, which 
we can write as P = P(T, Aq, A 2 ^ ). In other words, gauging the symmetry allows us to trade the 
dependence of pressure on the condensate phase for the dependence of pressure on the background 
gauge field. In a general gauge, P = P(T, —DQip, (D^tp) 2 ). Denoting \ = ^(D^ip) 2 , we have the 
pressure as a function of three thermodynamic variables: P = P(T, /u,x)- From it one can define 
the conjugate variables, 

dP = sdT + ndfi- f 2 d X ■ (9) 

Notice that s, n, and f 2 are functions of T, ju, and \. If the third term were absent from the Eq. ([9]), 
then s would be the entropy density, and n the charge density. We will give the interpretation of 
s, n and f 2 in the next section. 



III. HYDRODYNAMICS 



Let us set the external gauge field A^ to 0. The equations of ideal relativistic superfluid hydrody- 
namics are known. We will use the form written in In this formulation, the degrees of freedom 
are T, fj,, <p, and a unit four- vector satisfying •q^u^u 1 ' = — 1 (recall that we use the mostly plus 
convention for the flat metric tensor 77^). We will later identify with the velocity of the normal 
component. The equations consist of the conservations of the energy-momentum tensor T^ v and 
of the U(l) symmetry current j^, 



8^ = 0, 



(10) 
(11) 
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and a "Josephson equation" describing time evolution of (p, 

u^cp + [i = . (12) 

The stress-energy tensor and the current are expressed in terms of the hydrodynamic variables 
through 

= (e + P)u^ + Pff + f 2 d u <p , (13) 

f = nu» + fd^ip • (14) 

where the energy density e is defined by e + P = Ts + n/j,. Thus we have (d+2) hydrodynamic 
equations ([TO]) . (fTTj) . (fT2"j) for (d+2) variables T, </> and u^. One can derive from Eqs. (|TU|) — f|T^ 



d^su") = 0. (15) 

One can interpret this equation as the equation of entropy conservation (entropy is conserved since 
we are working at the level of ideal, nonviscous hydrodynamics). Thus s is the entropy density 
and is the velocity of entropy flow. In the two-fluid model only the normal component carries 
entropy; therefore u^ 1 is interpreted as the normal velocity. The two contributions to the current 
in Eq. (|14p can be interpreted as the normal and superfluid currents. Therefore n is the normal 
density, and f 2 is the analogue of the pion decay constant. (The equivalent of the superfluid density 
would be f 2 H-) 

Let us look at small fluctuations about an equilibrium state at fixed temperature, chemical 
potential and zero normal and superfluid velocities, i.e. we write T = Tq + T", fi = }1q + //, 

= (l,v 1 '), £j = diip = ££, where T', //, v i , and £' are small. In terms of £ the variation of 
pressure is dP = sdT+ (n+fif 2 )d^i - / 2 £d£, where £ = |£| . Further, let us assume that pressure is 
a smooth function of £ 2 at small £, so that <9P/<9£ = in equilibrium. The linearized hydrodynamic 
equations become (omiting subscript "0" on equilibrium quantities) 

' hP ■ 8 t T' + ^ dtfi' + f 2 + nd i v' i = 0, (16a) 



"ma + r S) d < r + (>w + v + w ** + /V m = ' (16b) 

w dtv'i + f 2 n dt£ + 8 d l T' + (n+nf 2 ) dm' = , (16c) 
d t ^ + d ifJ ! = 0. (16d) 

In this system of equations, pressure is taken as P = P(T, /U,£), and w = e + P is the density of 
enthalpy. The first equation is the linearized current conservation equation d^j^ = 0. The second 
equation is the linearized energy conservation c^T^ = 0. The third equation is the linearized 
momentum conservation, d^T^ 1 = 0. Finally, the fourth equation says that \x and £j are not 
independent because \x = —dtp while £j = d^ip. 

An interesting feature of the linearized hydrodynamic equations (|16jl is that they admit prop- 
agating mode solutions even if one ignores the fluctuations of the energy-momentum tensor, i.e. 
if one ignores the (normal) sound fluctuations. If one were to ignore the condition d fJ T^ v = 0, 
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together with fluctuations of temperature T' and velocity of the normal component v[, then the 
system (fT6|) becomes 

d 2 p 



dp? 

d t (l + d^' = o 

Fourier transforming all variables as e - lujt + lk - x ^ we f| nc [ a propagating mode with frequency 

/ 2 



(17a) 
(17b) 



a; 



2; 2 



a 2 p 



d 2 P' 

d 2 p' 
dfi 2 



> 0. 



(18) 



When expressing / 2 in terms of (d 2 P/d£ 2 ), we have assumed that P = 0(£ 2 ) at small £. The 
thermodynamic derivatives in Eq. (fTBj) are to be evaluated at £ = 0. The propagating mode (fTHJ) is 
the leftover of the second sound in superfluids which survives even if one ignores the daT^ u = part 
of the hydrodynamic equations. One expects that this mode is captured by the dual gravitational 
description which ignores the backreaction of the gauge fields on the metric. 

It is not difficult to find the eigenmodes of the full system (|16p . Taking all variables proportional 
to e ~ lu)t + lk ' x ; one finds a fourth order equation for frequency, 



bk 2 u 2 + ck q 







(19) 



where the coefficients a, b and c are independent of k. When b 2 > 4ac, the equation for uj 2 has 



v\k 2 . The first solution is the normal sound, and the 



two real positive roots oj 2 = v 2 k 2 and uj 2 
second solution is the second sound. In terms of thermodynamic derivatives, the coefficients a, b 
and c are given by 

2' 



Tw 



f d 2 p 
\dr 2 

f 2 Ts 2 



d 2 p \ 
or 2 ) 



& 2 p 

dp 2 



T(n 2 +wf 2 ) + 



/ d 2 p v 

\dTdfi) 



8 2 p 

dp 2 



Ts 2 



( 8 2 p 

\8Tdfi 



Ts n , 



(20) 

(21) 
(22) 



In particular, a, b, and c are positive. Let us now look at simple examples. In the non-superfluid 
phase f 2 = 0, and therefore v 2 = b/a, while the second sound is absent, v 2 = 0. Even in the 
non-superfluid phase, the speed of the normal sound looks like a complicated expression in terms 
of the derivatives of pressure P(T, fi). The expression for v 2 simplifies if instead of P(T, fi) we work 
with P(s,n). Indeed, the coefficient a is proportional to the Jacobian of the transformation from 
the (s,n) variables to the (T, p) variables. In terms of P(s,n), the speed of sound in the normal 
phase becomes 



n 
w 



dP 

dn 



s 

+ - 

w 



dP 

ds 



(23) 



Now that we have v 2 expressed in terms of P(s,n) without reference to the chemical potential, 
we can evaluate the thermodynamic derivatives in the canonical ensemble instead of the grand 
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canonical. In the canonical ensemble, the total charge (number of particles) N = nV is fixed, and 
therefore dn/n = —dV/V. It will be convenient to go from the (s,n) to the (S, e) variables. For 
the total entropy S we have dS /V = d{sV)/V = ds — s dn/n , while on the other hand the relation 
TdS = dE + PdV gives T dS/V = de — w dn/n . This implies that the speed of sound in the normal 
phase (f23j) can be written as 



\ de )s,N 

As another example, consider a scale-invariant theory, in which case pressure has the form P(T, u) = 
T d g(T/ fi), where g(T//j.) is a dimensionless function. The speed of normal sound evaluated from 
Eq. (fT9|) is v 2 = l/(d— 1), in either normal or superfluid phase. The speed of the second sound 
evaluated from Eq. (|19p can be written as 



vi = r 



1 + Ts 



un\ (d 2 P n + uf 2 d 2 P \ 



du 2 s dTdfi ) 



(25) 



We can see that in scale invariant theories, the cartoon expression (|18p can be recovered in the 
formal "large-entropy" limit Ts 3> un, together with s(§^?) 3> (n+/i/ 2 )( ^J^ ). In theories which 
do not have scale invariance, the speed of the second sound does not have the simple form (I25p . 
but can be easily determined from Eq. (|19p given the equation of state P(T,n). 

Note that the simple expression (|18[) ceases to be a good approximation to the speed of the 
second sound at T <C u. This is because at small temperatures we expect that the system can be 
described as a gas of massless Goldstone bosons, with entropy density s ~ T rf_1 . On the other 
hand, we expect in the low-temperature region that n ~ [i d ~ x which imples that the condition of 
large entropy can not be satisfied at T <C (jl. At low temperatures (high densities) the speed of the 
second sound has to be determined from the full equation (fT9|) , 



IV. DUAL GRAVITY DESCRIPTION 



On the gravity side, we study the Einstein-Maxwell theory with a complex scalar field. Fluctu- 
ations of the bulk metric correspond to fluctuations of T^ v on the boundary, while fluctuations of 
the f7(l) gauge field correspond to fluctuations of on the boundary. The scalar field is charged 
under the bulk U(l), and its background value corresponds to the condensate on the boundary. The 
mass of the scalar is a free parameter. We also need to specify the boundary conditions. The metric 
will be asymptotically AdS because we want to study the superfluid system using the AdS/CFT 
correspondence. The boundary conditions for the scalar correspond to the "normalizable" mode 
because we want to describe the system in which a charged operator has a vev. However, the 
boundary conditions for the U(l) gauge field should correspond to the "non-normalizable" mode 
because we want to describe the system in the background gauge field, as discussed above. In other 
words, we fix the value of at the asymptotic AdS infinity. Fixing the value of at infinity 
amounts to fixing the chemical potential in the boundary theory, and leads to a charged black hole 
in AdS. Fixing the value of A{ at infinity does not introduce additional "hair" for the black hole 
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because we expect the new black hole solution to be only metastable, in accord with field theory 
expectations. So we want to find a (metastable) stationary solution with the above boundary con- 
ditions, and then study small fluctuations around this solution, corresponding to hydrodynamic 
fluctuations in the boundary theory. 

In the following, we will ignore the backreaction of the gauge and scalar fields on the metric. It 
would be nice to include the backreaction, at least perturbatively to leading order. The action is 

1 

4? i M N ' 



S 



'd+l„ 



-9 



TP Z?MN 



+ {D u 4>){D M cj>f + m 2 



(26) 



where D M (f> = d M <p-iA M <p, F MN 
equations of motion are 



d M A N — d N A M , and capital Latin indices run from to d. The 



-9 9 



(27) 



-9 



MA„NB Tp \ _ „2 NA J 

f ab ) — e g J A 



-9 9 9 



(28) 



where the current is J A = i[(j>* (D A 4>) — <ft(D A <j))*\- We write the bulk scalar as <fi = --j^pe n and 
make a gauge transformation A M — ► A M + d M cp. In the new gauge, the phase <p disappears from 
the equations of motion, and the current becomes 



Ja 



P 2 A h 



(29) 



From the Maxwell equations (|28|) one can see that the background value for p would induce a 
(position-dependent) mass for the gauge field. This is the Higgs mechanism in the bulk. We take 
the (d+1) dimensional background metric to be of the following form: 



ds 2 = -yl 



-f( z )dt 2 +dx 2 + ^). 



(30) 



The metric with f(z) = 1 corresponds to pure AdS and f(z) = 1 — (z/z} l ) d corresponds to our 
black hole solution. To find the background solution for p and A M , we take all fields independent 
of t and x. The z-component of the Maxwell equations (j28[) now gives p 2 A z = which means we 
can choose A z = 0. The equations of motion become 



„d-l 



Z 



f 



yd-X 



P 



Aj 



m 



p- 



(31) 



yd— 3 



A{ 



e 2 

— f p 2 A t 

z 2 f 



(32) 



v d-3 



yd-3 



A' 



e 2 

-3 p 2 Ai . 
z z 



(33) 



In the limit when A% = they reduce to the coupled (At, p) system of equations studied recently 
in [3]. Note that we have a coupled system of non-linear ODEs — it would be interesting to see 
if the solutions may exhibit chaotic behavior. 



9 



The free energy 



Up to boundary counter terms, the free energy of the field theory is determined by the value 
of the action (|26|) evaluated on-shell, Q = —TS os + . . ., where the ellipsis denotes boundary terms 
that we will presently introduce. Employing the equations of motion, we may rewrite (|26l) as 



So 



d d x 



-99 



1 

+ 2 



■-h 



dz^jA^p 2 



(34) 



We have included a cut-off e because this on-shell action is naively divergent and needs to be 
regularized through the addition of boundary terms. It is difficult to treat the general case both 
succinctly and clearly and several full treatments already exist in the literature (see for example 



[la]). We proceed to regularize this action in the simple case d 
explicit form of the metric (|30p . the on-shell action reduces to 



3 and m 2 = —2. Using the 



So 



d d x 



tV ^A.A' v + 



tPP 



+ 



dz 



-gA^p 2 



The near boundary behavior of the fields takes the form 

A„ = (a^ + 0(z 2 )) + z{b^ + 0{z 2 )) , 
p = z(a + 0{z 2 )) + z 2 {b + 0{z 2 )) . 



(35) 



(36) 
(37) 



These boundary values have various reinterpretations in the field theory. For the gauge field, ao = p 
is the chemical potential while ai = — & is a superfluid velocity. Then bo ~ — n is proportional to 
the charge density while 6j ~ Jj are charge currents. For the scalar field, there exists an ambiguity 
171 ]. For a scalar operator 0\ of conformal dimension one, a = (0\) while b is a source. For a 
scalar O2 of conformal dimension two b = (O2) while a is a source. 

In regulating S os , we must carefully formulate the boundary conditions. For example, do we 
wish to keep A^ or A'^ fixed on the boundary? Keeping Aq fixed corresponds to keeping the 
chemical potential fixed and thus working in the grand canonical ensemble in the field theory. 
In the spatially homogenous case where we can set A z = 0, varying Am in the bulk leads to a 
boundary term proportional to A'^SA^. Without an additional boundary term, we are working in 
an ensemble where a M is held fixed. If we would like to work in the canonical ensemble, at fixed 
charge, it is b$ that must be held fixed at the boundary (5A' = 0). To accommodate this change, 
we would need to make what amounts to a Legendre transform and add the boundary term 



i J d d xA A' 



-pQ/T 



(38) 



to the action (|26p . We will work in an ensemble where is held fixed and thus need no such 
further boundary terms. 

A similar decision needs to be made about the scalar operator. It is most natural to work in 
an ensemble where the value of (Oi) is fixed instead of the source for the operator. However, we 
still must decide whether we want our scalar operator to have conformal dimension one or two. As 
the ensemble where the source for 0\ is fixed is equivalent to the ensemble where (O2) is fixed, we 
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enow that the ensembles where (Oi) and (O2) are fixed must be related by a Legendre transform 



12|. 



Consider first the case where (Oi) is fixed. The on-shell action is naively divergent, and we 
must add a counter-term. The counter-term and the pp' term in the on-shell action combine to 
give 



1 



1 



2^ PP ~ 2^ P 



\ab + 0(e) . 



(39) 



In the Legendre transformed case, where (0 2 ) is fixed, we add two counter-terms, one to ensure 
that we hold bp' fixed at the boundary instead of bp and one to control the divergence: 



1 



1 



1 



2? PP ~ * PP + 2^ P 



\ab + 0(e) . 



(40) 



Recalling that this regularized on-shell action is the negative of the free energy and assuming a 
spatially homogenous system so that we may divide out by a factor of the volume V, we find for 
the free energies that 



ni(p,,(,Oi)/v 



-pn + i- J +(-l) l O l 2 - A il A u g^p 2 ^-gdz 

Jo 



(41) 



Recall Q = —PV. Having fixed Oi, p and £, the conjugate boundary values CijOj, Q, and J s are 
then determined through the dynamics of the gravitational theory. 

Since Oi(Oi) and ^2(^2) = ^i(Oi) + O1O2V are Legendre transforms of each other, we have 



1 _ 1 dn 2 

V oO\ V 0O2 



(42) 



Thus critical points of the free energies correspond to gravitational solutions where at least one of 
the two Oi vanish. 



Numerical results 

The nonlinear differential equations (j31H33[) appear to be intractable analytically. However, 
it is relatively straightforward to integrate the equations numerically. The results of this section 
are most succinctly summarized by the two phase diagrams for the scalars 0\ and O2 shown in 
Figure [H 1 

To see where these phase diagrams come from, we begin by reviewing the case £ = studied 
in the canonical ensemble in [3] while here we choose to work at fixed p. Given £ = 0, the third 
differential equation (f33|) drops out. Plots of the expectation value (Oi) versus temperature are 
shown in Figure [21 as the black curves on the far right. At high temperature, (Oi) = 0, but at 
the critical temperature T c , there is a second order phase transition where the expectation values 



1 We set e = 1 in this section. 
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a) Mb) M 

FIG. 1: The phase diagrams for the theory with a scalar with a) conformal dimension one and b) conformal 
dimension two. The solid blue line indicates a second order phase transition while the solid red line (in 
between the dashed lines) indicates a first order phase transition. The dashed blue lines are spinodal curves, 
while the red dot indicates the tricritical point. 

become nonzero. For Oi, T c = 0.213/x while for O2, T c = 0.0587/z. Near but slightly below T c the 
scalars exhibit the standard mean field scaling with the reduced temperature 

(Oi) ~ (T c — T) 1 / 2 . (43) 

The most straightforward way to see that the phase transition is second order is to examine a plot 
of the free energy versus temperature: Qi is smooth at T c . Using eq. ([4"Tj) . we have produced Figure 
[3^ for the scalar 0±, which indeed shows this smooth behavior. We do not show a very similar plot 
for the second scalar 02- 




a) V^ 2 -f 2 b) V^-^ 2 

FIG. 2: The condensate as a function of temperature for the two operators: (a) 0\ and (b) Oi- The curves 
in the plots, from right to left, are for = 0, 1/4, 1/3, 2/5, and 1/2. 

A more elaborate demonstration that the phase transition is second order comes from an in- 
vestigation of Qi as a function of the order parameter Oi near T c , recovering completely standard 
results in the Landau-Ginzburg mean field theory of phase transitions. We have found numerically 
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FIG. 3: The difference in free energy Afij between the phase with a scalar condensate and without one as 
a function of T/fj.; a) £ = and b) = 4/7. 

that 

= i( HwJ +A( } lwJ (44) 

fits the free energy curves extremely well. Moreover, on and are nearly linear in T: 

a t = 3.07(T - T c )/T c , fa = 0.743 - 0.899(T - T c )/T c , (45) 
a 2 = 5.13(T - T c )/T c , f3 2 = 1.34 - 1.45(T - T c )/T c . (46) 

By definition, a, passes through zero at the phase transition. 

As we increase the superfluid velocity £, nothing dramatic happens immediately. The phase 
transition remains second order although T c decreases as can be seen from Figure [2j Because 
the decrease is due to the additional kinetic energy of the system, we expect the decrease to be 
quadratic in £, which is born out by the shape of the second order lines in the phase diagram, 
Figured) Numerically, we find that T c (£) T c (0) - A£ 2 //" where A « 0.27 for O x and A » 0.14 for 

o 2 . 

However, there exists a critical £ above which the phase transition becomes first order. For 0\, 
this critical velocity is £ = 0.342^ while for Oi it is £ = 0.274^. From Figure it is clear that 
something interesting must happen because the curves (Oj) become multi- valued for sufficiently 
large £. It is possible to see that the phase transition is first order in different ways. The simplest 
is to look at the free energy as a function of temperature. Figure Eb presents this classic swallow 
tail shape for 0\ and = 4/7. At the phase transition, the free energy is continuous but not 
differentiable. The two nonanalytic points in the free energy curve are "spinodal points" or points 
beyond which one of the phases ceases to exist even as a metastable minimum of the free energy. 

To demonstrate more convincingly that the phase transition becomes first order, we computed 
the free energy as a function of the order parameter near the putative tricritical point. We found 
numerically that the free energy is well described by the sixth order polynomial 
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1.06 



T T 

a) Mb) P 

FIG. 4: The speed of second sound as a function of T / /z, computed by evaluating thermodynamic derivatives 
in Eq. (fT8|) : a) 0\ scalar, b) O2 scalar. The speed of second sound vanishes as T — > T c and appears to 
approach a constant value as T — ► 0. 



At the tricritical point, a, and both vanish. Moreover, near the tricritical point, on and vary 
linearly with T and £ while 7, is positive and roughly constant. 

Before moving to a discussion of second sound, we would like to point out one nice feature of 
Fig. [2j The curves approach each other at low temperature. This agreement is a consequence of 
the fact that at T = 0, the pressure P can only be a function of the Lorentz invariant quantity 



2 



£ 2 , as discussed in Section HH 



The speed of second sound as a function of temperature is shown in Fig. [H The plots were 
computed in the probe approximation using Eq. (|18|) . We have set the superfluid velocity £ = 0. 
The behavior close to T c is qualitatively similar to that of superfluid 4 He 18]. As the temperature 
is decreased from T c , the speed of second sound rises rapidly from zero and eventually levels off. 
Experimentally, it is difficult to go to very low temperature and remain within the hydrodynamic 
limit. The scattering length for the phonons approaches the system size. Numerically, we have 
also had difficulty finding solutions at very low temperatures; the curves in Fig. U] terminate where 
our numerics fail. Theoretically, however, the exp ectation for 4 He is that second sound at T ~ is 
a sound wave supported by a gas of phonons [19]. Thus, the speed of second sound close to T = 
should approach v\ = v 2 s /(d — 1). 

The naive extrapolation of the curves in Fig. H] suggest that v\ is 1/3 in the 0\ case and 1/2 in 
the O2 case. The fact that the curves in Fig. [J] level off at low temperature is caused numerically 
by a similar growth in the susceptibility and the pion decay constant in Eq. (|18p . That we work in 
the probe approximation suggests our model may not be reliable at low temperature anyway. We 
have neglected the coupling between our Abelian-Higgs sector and the metric. The Abelian-Higgs 
sector on its own does not support ordinary sound. 
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V. CONCLUDING REMARKS 



We have discussed the hydrodynamics of relativistic superfluids. Moreover, we presented a 
holographic model that reproduces many of the familiar features of superfluids, including a critical 
superfluid velocity above which the system returns to its normal phase and a second order phase 
transition at zero superfluid velocity between the normal and superfluid phases. 

In many respects, our holographic model is similar to the Landau-Ginzburg mean field treatment 
of phase transitions. Near the phase transition, we saw that the potential function for the order 
parameter has the familiar polynomial expansion V(<p) = acj) 2 + (3^ + . . .. In our case, however, the 
values of a and (3 and indeed the entire structure of V(cp) are not set directly by the field theorist 
but are instead encoded in a nontrivial fashion by the bulk gravitational solution. By assuming a 
simple gravitational model, we are led to a particular V(4>) that is valid not just near T c but at all 
temperatures in the superfluid phase. Moreover, from the model one can extract not only the static 
thermodynamic properties, but also quantities relevant for time-dependent, dynamic processes, like 
the kinetic coefficients and the correlation functions. 

In the future, it would be very interesting to extend the numerical results above. Two obvious 
directions present themselves. In the calculation of the speed of second sound from derivatives of 
the pressure P, one could go beyond the probe approximation and include the back reaction of 
the Abelian-Higgs sector on the metric. We hope that the resulting speed of second sound will be 
phenomenologically meaningful not only near T = T c but all the way down to T = 0. 

Another interesting project for the future would be to find the sound wave poles in the density- 
density correlation function. By studying fluctuations of the scalar and gauge field of the form 
e —iuit+ikx j n fog probe limit, one should be able to isolate a pole of the form 1/(lo 2 — v\k 2 + . . .) 
in the Fourier transform of the retarded Green's function for the charge density. Here the ellipsis 
denotes higher order (damping) terms in k. Beyond the probe limit, there should also be a pole 
corresponding to the propagation of ordinary sound. 



Note Added — While we were completing this paper, we learned of ref. 26|] which has some overlap 
with this work. 
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APPENDIX A: COMPARISION WITH THE CARTER-KHALATNIKOV-LEBEDEV 
FORMULATION OF SUPERFLUID HYDRODYNAMICS 



In this Appendix we show that the hydrodynamic equations written in Sec. IIIII are equivalent 



to the set of equations proposed previously by Israel, Carter, Khalatnikov, and Lebedev [2CJ, 



Refs. 



21 



22j, |23 . \24 , |25H . The formulations of Ref. 20| and Ref . 21[ have been shown to be equivalent in 



23H. We will follow the notation of Ref. [2c 
In Ref. [23], superfluid hydrodynamics is formulated as follows. First, one postulates that the 
thermodynamic properties of the superfluid are defined by a scalar function A, which is a function 
of two 4- vectors s" and n", which are the entropy density and the particle number density. Since 
A is a Lorentz scalar, the number of variables that it depends on is three: s^s^, n^n^, and s M n M . 
One defines two Lorentz vectors and fi" from 

dA = Q^ds" + [i^dn" 

The superfluid hydrodynamic equations are 

d„.s 



(Al) 



TV 



, (A2) 
: , (A3) 
d^Hu ~ dvHn = , (A4) 

s"(0 M e„ - a v ej = . (A5) 

At the first sight, Eqs. (IA2h - ()A5j) do not bear much resemblance to Eqs. (|10p - (|15p . However, 
these systems of equation are in fact equivalent. The variables appearing in the Carter-Khalatnikov 
formulation can be identified with the hydrodynamic variables used in this paper as in the following 
relations, 

A = e-f 2 (d^) 2 , (A6) 
s" = su" , (A7) 
nu" + fd^cp , (A8) 
-dp<p , (A9) 

(A10) 

We now show that any solution to the hydrodynamic equations ()LQ[) — (jL5|) satisfies Eqs. (|Aip - (|A5p . 
upon the subsitutions (|A7^ - (fA"l0j) . First, Eqs. JA2]) and (|X3]) coincide with Eqs. JTS]) and (JIT]) due 
to Eqs. ([ATj) and (JM]). Furthermore, Eq. (fA"4|) is trivially satisfied by Eq. (JA~9]) . 
Now let us check Eq. (|Aip . The left hand side is 

dA = de- d(f(d fl ip) 2 ) = Tds + ndn - d{f 2 {d"^)) , (All) 

where the thermodynamic relation e + P = sT + \xn and Eq. have been used, while the right 
hand side is 

- [-(Ts + finjup + nd^Lp] (u^ds + sdu") - d^[u"dn + ndu" + d(fd^ip)] 

Et} dn + -i^B.m ds - n"8..cn dn - B. .in d( f 2 8"cn\ (A12) 



Mm 



- [-{Ts + ^n)u^ + nd^ip] 



7 + ds + -u^d^cp ds - u^d^ip dn - d„(p d(fd^ip) 
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where we have used u 2 = —1, u^du 11 = 0. Now from the Josephson equation u^d^ip = —\i one sees 
immediately that (|A12|) is identical to (jAlip . Thus, Eq. (|A1|) is verified. 
The last equation that has to be checked is Eq. (IA5I) . We first write 



where = d^su^ 



has been used. We expand 



where Eq. (fTT|) has been used; 



fin 



n 



T H )ua + -<9«v? 



n 



sd v T + nd v u + sfid u - + su^d^ipdv - + nu^d^dyip 



n 



Combining Eqs. (IA14I) and (jA15|) . using the Josephson equation, one finds 

s»{d^Q u - 3 U Q^) = -d^[(Ts + nn)u^u v ) - <9 M (/ 2 cKV)cW - sd u T - nd uf i . 



(A13) 



(A14) 



(A15) 



(A16) 



It is easy to check that the right hand side is equal to d fl T fl u up to an overall sign, with T^ u defined 
in Eq. (|13p . and hence is equal to zero. 
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